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1.  Introduction 


The  electron  spin  in  a  single  quantum  dot  is  one  of  the  perspective  realizations  of  a  qubit  for  the 
implementation  of  a  quantum  computer.  During  the  last  decade  several  control  schemes  to 
perform  single  gate  operations  on  a  single  quantum  dot  spin  have  been  reported  (1^4).  Here  we 
propose  a  scheme  that  allows  performing  ultrafast  arbitrary  unitary  operations  on  a  single  qubit 
represented  by  the  electron  spin.  The  idea  of  geometric  manipulation  of  the  qubit  wave  function 
has  been  recently  developed  into  a  new  research  direction  called  geometric  quantum 
computing  (15-17).  The  main  motivation  of  this  development  is  the  robustness  of  geometric 
quantum  gates  against  noise  (15,  18,  19).  In  this  report,  we  demonstrate  how  to  use  the 
geometric  phase,  which  the  Bloch  vector  gains  along  the  cyclic  path,  to  prepare  an  arbitrary  state 
of  a  single  qubit.  We  show  that  the  geometrical  phase  is  fully  controllable  by  the  relative  phase 
between  the  external  fields.  Realistic  implementation  of  the  proposed  design  using  the  electron 
spin  in  a  charged  quantum  dot  as  an  example  of  a  qubit  is  discussed. 


2.  General  Equations  of  Motion 


Let  us  consider  the  coherent  Raman  excitation  in  the  three-level  A-type  system  consisting  of  the 
two  lowest  states  of  electron  spin  |0)  =  |  —  X)  and  |0)  =  \X)  coupled  through  an  intermediate 
trion  state  | T)  consisting  of  two  electrons  and  a  heavy  hole  (5)  (figure  1).  We  assume  that  the 
trion  state  is  far  off-resonance  with  the  external  fields  to  ensure  that  decoherence  on  the 
trion-qubit  transitions  can  be  neglected.  The  electron  spin  states  are  split  by  an  external  magnetic 
field;  the  separation  energy  is  hue.  The  total  wave  function  of  the  system 

\*(t))  =  aQ(t)\0)  +  a1(t)\l)  +  b(t)\T),  (1) 


where  a0,i(£)  and  b(t)  are  the  probability  amplitudes,  is  governed  by  the  time-dependent 
Schrodinger  equation  with  the  Hamiltonian 


H 


0 

0 

—  (Qp(t)  +  Cls(t)) 


0 

LUe 

—(Ctp(t)  +  f ls{t)) 


—(Qp(t)  +  £ls(t))  ^ 

—  (Clp(t)  +  Qs(t))  ) 

Up  ) 


(2) 


where  f Ip, sit)  —  Qpo,so(t)  cos[c up,st  +  (f>p,s{t)],  &p,s(t)  —  Qpo,so(t)  cos[uptst  +  (f>p,s(t)]. 
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ttPO,so(t)  =  tk}T.T\ EpS{t) / h,  and  fipo,so (t)  =  i-i\T,ToERS(t)/h  are  the  Rabi  frequencies,  / x0t,ti 
are  the  dipole  moments,  EP^s(t)  are  the  pulse  envelopes,  ujp.s  are  the  center  frequencies,  Op. sit) 
are  the  time-dependent  phases,  and  Hujt  is  the  energy  of  the  trion  state. 


Figure  1 .  Energy  structure  of  the 
three-level  system 
comprised  of  the  two 
electron  spin  states  and 
the  trion  state. 


We  are  considering  here  a  case  of  linearly  chirped  pulses  such  that 

=  (pp,s  +  «p,5f2/ 2 ,  (3) 

where  4>P^  are  the  initial  phases  and  ap.s  are  the  chirps  of  the  pulses. 

Using  transformation  | ']/(£))  =  UpwA^it)),  where 


/  1 


UrwA  — 


0 


Q  ei(us-uP)t 


0  \ 
0 


\0 


0 


D — iujpt 


(4) 


IT-1 

URWA 


(  1 


0 


0  \ 


0  e-i(ws-wp)t  Q 


(5) 


\  °  o 

we  make  the  rotating  wave  approximation  (RWA)  by  neglecting  the  rapidly  oscillating  terms  with 
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frequency  2c os,  2up  and  lus  +  oJp-  In  the  RWA,  the  Hamiltonian  has  the  following  form 


H  U pyy  j^HU RW  A 


ihuRWAuRWA 


flp,  +  fl%,eiAuJt 


n*p+e~ 


+ 


/  0  0 
+  h  o  — Acu 
Vo  o 


n*p+  +  Q.*s+el 


n*p+e 


2(Ac 0  —  LUe) 

e-iAut  + 


Hp+  +  Vts+c~lAult 
ClP+(t)elAuJt  +  VLS+ 
— 2c Op 


fip+  +  fis+e  lAuJt 
ClP+(t)elAuJt  +  f }s+ 
— 2AP 


where  Ap  —  c op  —  cup,  Acu  —  c up  —  lus,  flp+  — 
f2p+  =  f2p0(f)e*'?ip^V  0p+  =  05o(f)e*'?is'^V 

Assuming  large  detunings  of  the  pump  and  Stokes  field  frequencies  from  the  transition 
frequencies  to  the  trion  state,  we  apply  the  adiabatic  elimination  of  the  trion  state.  In  that  case, 
b  ~  0  and  we  find 


b  =  ^  [(sit+  +  siJA  )  a ,(()  +  (ncp-A"'  +  n%+)  a,(t)] . 


Substituting  equation  7  into  equations  for  «0: 1  after  some  algebra,  we  obtain  the  following  form  of 
the  Hamiltonian  for  the  effective  two-level  system 


jj  =  _h  (  ^aco(f)  \  8 

2  ^  n*eff(t)  nacl(t)-25  )  ’ 

where  <5  =  Lue  —  Auj, 

Maco(t)  =  — (Opo(f)  +  fiso(t)  +  2Hp0 (t)Qso(t)  COS  [<j>P(t)  -  4>s(t)  +  Ac vt]  ,  (9) 

ZL \  p 

ttacl(t)  =  7TT—  (^so(0  +  &Po(t)  +  2fiso(t)fip0(f)  COS  [(f) P(t)  ~  <t>S(t )  +  Awt]  ,  (10) 

ZL \  p 
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^—iAujt 

^ eff(t )  =  -— r (fipo(t)Opo(t)  +  figo (£)Og0 (t) 

ZL \  p 

+  nPO(t)nSo(t)ei[,t,p{t)-(l,s{t)+Au,t]  +  (W^so(*)e_#pW_^w+A^])  . 

Since  f2po(f)fiso(£)  =  Hpo(t)flso(£),  we  have 

^—iAcot 

^ eff(t )  =  ~PTT (fipo(t)Opo(t)  +  Ogo  (t')foso(t') 

ZLA  p 

+  2f)p0(f)f)go(£)  cos[0p(t)  -  +  Acot])  . 


(11) 


(12) 


Making  the  transformation  TU))  =  U\x\>(t)),  with 


U  = 


3— i8t./2+i 

0 


0 

j~i5t/2+i  dt'Clst(t') 


—  e~iSt/2+i  dt'nat(t')  j  _ 


where  Clst{t)  =  (f)ac0(£)  +  fiaci(f))/4,  we  can  rewrite  equation  8  in  more  symmetric  form.  In  the 
new  basis,  the  Hamiltonian  takes  the  form 


H  =  U~XHU  -  im~lU  =  H  +  h(nst(t)  -  5/2)1 
h  f  6 +  ndif(t)  neff(t) 

^e//W  ~  ^ dif{t ) 


(14) 


where  fidi/(t)  =  (nac0(f)  -  Oaci(t))/2. 

Taking  into  account  the  exact  time-dependence  of  the  phases  in  equation  3  we  apply  another 
transformation  |  (&(£))  =  Wpl'k(f)),  with 


Uf  — 


opt1  /  4  Q 

0  e-<*2/4 


=  Ss"* 


(15) 


where  (  =  aP  —  as-  Therefore,  in  the  field  interaction  representation  the  Hamiltonian  takes  the 
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following  form 


H  =  Up' HU F  -  itiUp'Ur 

_h  (  5  +  ndif(t)  vteff(t)e-^2u\  hf(t  0  \ 

2  ^  niff{t)e*w*  -5-ndif(t)  y  2  ^  0  -C*  J 

_h(  S-(t  +  f hifit)  fie//(f)e-<t2/4  \ 

2  y  n*eff(t)e^t2/4  -5  +  (t-ndif(t)  J 

=  _h(  S(t)  Cleff(t)eiA *  \ 

2  ^  ti*eff(t)e-iA*  -S(t)  J  ’ 


where  5(f)  =  5  +  (as  -  otp)t  +  Qdif(t),  A0  =  (j)P  -  (ps, 


^ eff(t ) 


topo(t)toso(t)  A  Qpo(t)Qpo(t)  +  Qso(t)&so(t) 

2A  p  \  Qpo(t)Qso(t) 

npo{t)nSo{t)  -2i[Acf>+Aujt+(t2 /2]\ 

Up0(t)nSo{t)  J 


e-i[A<f>+Aut+(t2 /2] 


(17) 


In  general,  the  differential  AC  Stark  shift,  Vtdif(t)  is  not  zero  and  has  to  be  taken  into  account. 
Using  the  definition  of  the  Rabi  frequency,  we  can  rewrite  equations  9  and  10  in  the  form 

naco(t)  =  ^x;(EUt)  +  E2s(t )  +  2 Ep(t)Es(t)  cos  [Mt)  -  Mt)  +  A ut]  ,  (18) 

toaci(t)  =  2 h?Ap(E2p^  +  +  2UP(f)Us(f)  cos  [<f>P(t)  -  Mt)  +  Awf]  .  (19) 

In  the  case,  when  //0t  ~  Hi t,  the  AC  Stark  shifts  are  the  same  for  both  qubit  states,  so  that 
fl dif(t )  =  0,  and  we  have 


H  =  --(  \ 

2  V  V-:„(t)e-‘A*  -6  +  0  ) 

Note  that  we  do  not  require  completely  overlapped  pulses  here. 


(20) 
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3.  Non-Impulsive  Case 


In  some  excitation  schemes,  due  to  selection  rules  that  take  into  consideration  the  polarization  of 
the  external  field,  the  pump  and  Stokes  field  interact  only  with  the  corresponding  transitions  and 
the  general  Hamiltonian  can  be  simplified.  For  that  case,  we  can  obtain  the  correct  form  of  the 
Hamiltonian  by  putting  flpo(t)  =  &so(t)  =  0  in  equations  9,  10  and  17.  It  results  in 


Qaco(t)  2,Ap  (^)  ’ 

(21) 

^acl(t)  ~  2^p^So(t)  > 

(22) 

(t)  ~  ^so(^)]  ) 

(23) 

(24) 

Therefore,  in  the  field  interaction  representation,  the  Hamiltonian  has  the  form 


ft  /  i(t)  neff(t)e‘A* 
2  l  n;„(t)e-‘A*  -S(t) 


where 


5(t)  =  5  +  (as-  aP)t  +  [Hp0(f)  -  ft|0(i)]/ 4A; 


4.  Adiabatic  Solution 

The  Hamiltonian  in  equation  25  controls  the  dynamics  of  the  qubit  wave  function  in  the 
approximation  of  the  adiabatic  elimination  of  the  trion  state.  Here  we  consider  the  adiabatic 
excitation  of  the  qubit  and  find  the  adiabatic  solution  of  the  Schrodinger  equation  with  the 
Hamiltonian  in  equation  25. 

Since  the  phase  factor,  e',  A':>,  of  the  coupling  term  in  equation  25  is  time  independent,  it  is 
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convenient  to  use  the  following  transformation  |<f>(f))  =  A\'-V(t  j),  where 


A  =  |0)(0|  +  eiM\l)(l\  =  eiM/2e-iA^z/2  ? 

so  that  the  new  wave  function  is  governed  by  the  Hamiltonian 

H  =  AHA  1  =  --(  6{t)  °e(t)  ) 

2  ^  ne(t)  -5(t)  J 

h 

=  --  (S(t)crz  +  Qe(t)crx)  . 


(27) 


(28) 


To  solve  the  Schrodinger  equation  in  the  adiabatic  representation,  we  apply  another 
transformation:  | $(£))  =  R(t)\&(t)),  where 


m 


cos  6{t)  sin  6{t) 
—  sin  9{t)  cos  0(f) 


—  eiS(t)cry 


(29) 


and  tan 1 20 (7;)]  =  Vte(t)/5{t).  In  the  new  basis  the  Hamiltonian,  equation  28,  takes  the  form 


H(t) 


R(t)H(t)R~'(t) 

-HX(t)crz/2, 


0  ) 

fv^w+nfw  j 


(30) 


where  A (f)  =  \/82{t)  +  f^(f). 

As  we  see,  the  Hamiltonian  in  equation  30  is  diagonal  in  the  adiabatic  basis  and  we  can  readily 
write  down  the  solution.  However,  since  the  transformation  R(t )  is  time  dependent,  an 
additional  nonadiabatic  coupling  term  is  present  in  the  general  Schrodinger  equation 


ih\$(t))  =  ~h\(t)crz\$(t))  -  M(t)cry\$(t)) , 


(31) 


where 


Qe(t)5(t)  -  S(t)Qe(t) 

2  MW+m) 


(32) 


Neglecting  the  nonadiabatic  coupling  term  in  equation  31,  we  readily  obtain  for  the  qubit  wave 
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function  in  the  original  basis 


|  'T(f))  =  e^c 


0-id(i)CTy  pi^-CTz 


,i0(O)<Ty  p-i^r<Tz 


1^(0)). 


(33) 


where  Kit)  =  J'j  dt'X(t').  Note,  that  the  general  form  of  the  evolution  operator  in  equation  33  is 
well  justified  if  the  following  condition  |  Qe(t)S(t)  —  Vle(t)5(t)\  <C  A  3(f)  is  valid. 

In  the  case  of  completely  overlapped  pulses,  fiP0(t)  =  Qso(t),  with  identical  chirp  rates, 
ap  =  as,  for  the  resonant  qubit,  8(t)  —  6  —  0,  we  have  9(t)  =  9(0)  =  7t/4  and  the 
transformation  matrix  becomes  R(t)  =  R( 0)  =  e*1'7^4.  Therefore,  the  unitary  evolution 
operator  for  the  wave  function  of  the  resonant  qubit  takes  the  form 

=  /  cos  (S(t)/2)  ieiA(t>  sin  (S(t) / 2)  \ 

Y  ie_*A<^sin  (S(t)/2)  cos(S(t)/2)  J 

=  cos  (S(t)/ 2)  I  +  isin  (S(t)/2)  (elA<V+  +  e~lA<t>(T^) 

=  cos  (S(t)/2)  I  -  i  sin  (S(t)/2)  (n  ■  a)  =  e~iS{t)n  (Tl2  ,  (34) 


where  S(t )  =  /J  dt'Qe(t')  is  the  effective  pulse  area,  n  =  (—  cos  A  (j),  sin  A  0, 0), 
cr±  =  (crx  ±  icTy)/ 2  are  the  Pauli  raising  and  lowering  operators,  and  <T:r  .y.z  are  the  Pauli 
operators.  Note,  that  the  nonadiabatic  coupling  term,  equation  32,  is  zero  for  the  resonant  qubit 
and  the  solution  of  the  Schrodinger  equation  in  the  adiabatic  approximation,  equation  34,  is  the 
exact  solution. 


The  density  plots  of  the  population  and  coherence  (|a0(T)aj(T)|)  at  the  final  time  (after  the  pulse 
excitation)  as  a  function  of  the  effective  pulse  area,  S(T),  and  the  dimensionless  frequency  chirp 
parameter,  a' /tq  ,  described  by  the  unitary  evolution  operator  in  equation  34  is  depicted  in 
figure  2.  We  use  the  Gaussian  shape  for  the  pulse  envelopes,  assuming  that  linear  chirp  is 
obtained  by  applying  a  linear  optics  technique,  meaning  that  a  transform-limited  pulse  of  duration 
r0  is  chirped,  conserving  the  energy  of  the  pulse  (6,  7).  The  temporal  (a)  and  spectral  (a’)  chirps 
are  related  as  a  =  a/r0_4/(l  +  a/2/r(f  )  (6,  7),  where  r0  is  the  transform-limited  pulse  duration. 

We  observe  the  Rabi  oscillation  regime,  when  the  population  of  the  qubit  states  is  changing 
between  0  and  1  while  the  coherence  is  changing  between  0  and  1/2.  This  behavior  does  not 
depend  on  the  chirp  rate,  since  the  effective  Rabi  frequency  Qe(t)  is  determined  by  the  product  of 
the  pump  and  Stokes  Rabi  frequencies:  flpo,so(t)  =  o  exp {— t2/ (2r2)}/[l  +  a'2 /tq]1/4  with 
the  chirp-dependent  pulse  duration  r  =  r0[  1  +  a:/2/r(j] l,/2  and  the  amplitude  (6,  7). 

In  turn,  for  the  off-resonant  qubits,  <5^0,  the  evolution  operator  in  the  adiabatic  approximation 
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Figure  2.  The  density  plot  of  the  1 1)  state  population  (a)  and  coherence  (b)  as  a  function  of  the 

effective  pulse  area  and  frequency  chirp;  ap  =  as,  5  =  0.  Initially,  only  the  |1)  state  is 
populated. 


takes  the  form 

U(t ) 


where 


e*£(t)  cos 

sjn 


<F)e*A0si \ 

cos  0(t)  J 


cos  9{t) 


sin  9{t ) 


1 

71 


i 

71 


5 

W+W) 

6 

7^+W) 


? 


and  £(f)  =  \  f*  yj52  +  f 11(t')dt'  is  the  effective  pulse  area. 

Note  that  for  the  off-resonant  case  the  transformation  matrixes  are 


R~\t) 

R(  o) 


—  sin  9(t) 
cos  9  [t] 


1 


(35) 


(36a) 

(36b) 


(37) 


Figure  3  demonstrates  excitation  of  the  off-resonant  qubit.  As  expected,  the  population  of  the 
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(a)  (b) 
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Figure  3.  The  density  plot  of  the  1 1)  state  population  (a)  and  coherence  (b)  as  a  function  of  the 
effective  pulse  area  and  frequency  chirp;  ap  =  as,  Sto  =  0.75.  Initially,  only  the  |1) 
state  is  populated. 


off-resonant  qubit  at  the  final  time  is  not  changed  by  the  external  fields  as  long  as  the  pulse 
excitation  parameters  are  in  the  adiabatic  regime.  This  is  the  regime  of  adiabatic  return. 
However,  we  still  observe  the  Rabi  oscillation  for  the  value  of  the  chirp  \a'\  <  5tq  .  This  is  the 
area  of  nonadiabatic  population  transfer  where  the  nonadiabatic  coupling  term  cannot  be 
neglected  in  equation  3 1 . 


5.  Bloch  Vector  Representation 


The  dynamics  of  the  qubit  wave  function  can  be  described  equally  well  using  the  Bloch  vector 
representation.  In  addition,  the  Bloch  vector  formalism  allows  a  very  nice  and  clear  geometrical 
interpretation  of  qubit  dynamics  (8).  In  this  section,  we  give  a  short  overview  of  the  Bloch 
picture. 

A  general  state  of  a  qubit  can  be  described  as 

|'T)  =  ao(t)|0)  +  ai(t)|l)  =  cos  (P/2)  |0)  +  e*asin(^/2)  |1) ,  (38) 

where  a  and  (3  are  the  phase  parameters.  Up  to  an  insignificant  global  phase,  the  wave  function 
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can  be  mapped  into  a  unitary  Bloch  vector  B  =  (u,  v,  w),  as  shown  in  figure  4. 


Figure  4.  The  Bloch  vector  representation  of  the  qubit  state. 

Excitation  of  the  qubit  by  an  external  field  corresponds  to 
the  rotation  of  the  B  vector  about  the  pseudo  field  vector, 
12.  with  components  determined  by  the  effective  Rabi 
frequency  f2e(i),  detuning  5,  and  the  relative  phase  Ac j>. 


To  use  Bloch  vector  representation,  we  construct  the  qubit  density  matrix 


Q  =  l*><*l  = 


000  001 
010  011 


1  /  1  +  cos  f3  e  ia  sin  f3 

2  l  eia  sin  f3  1  —  cos  /3 


(39) 


where  Qij  =  ai(t)a*(t),  i,j  =  0, 1.  Using  a  Pauli  matrix  decomposition 


Q  =  -( I  +  B  a ) 


(I  +  U  ■  (Tx  +  V  ■  <Ty  +  W  ■  (Tz) 


=  —  (T  T  cos  a  sin  /3  ■  crx  +  sin  a  sin  (5  ■  <ry  +  cos  f3  ■  crz)  , 


(40) 


we  identify  relation  between  the  components  of  the  qubit  wave  function,  the  qubit  density  matrix 
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elements,  and  the  Bloch  vector  components.  Therefore,  we  obtain  the  following  expression 


( " ) 

(  £>01  +  010  ^ 

/  cos  a  sin  / 3  \ 

B  = 

« 

= 

*(001  —  0io) 

= 

sin  a  sin  (3 

{w  ) 

R5 

o 

o 

1 

\  cos  /3  ) 

(41) 


Taking  into  account  the  Hamiltonian  in  equation  25,  the  equation  of  motion  for  the  density  matrix 
ihg  =  H q  —  qH  takes  the  following  explicit  form 


*000  = 

(e  iA</,0oi  -  eiA</,0iO)  Ue(f)/2  , 

(42a) 

*001  = 

— A (t) 0oi  —  (011  —  000)  Ue(f)e*A<^/2  , 

(42b) 

*010  = 

S(t)gio  +  (011  —  000)  Ue(t)e  tA<^ / 2, 

(42c) 

*011  = 

-  (e“iAV  -  eiA^0i 0)  ne{t)/ 2. 

(42d) 

Using  the  relations  in  equation  41,  we  can  also  write  the  dynamic  equation  in  the  Bloch  vector 
representation  as 


§  (■ t )  +  wQe  (t)  sin  A (f) , 

(43a) 

- u5(t )  +  w£le(t)  cos  A 0 , 

(43b) 

-uQe  (t)  sin  A0  —  vQe  (t )  cos  Acft . 

(43c) 

Introducing  a  pseudo  field  vector,  ft,  with  components  determined  by  the  effective  Rabi 
frequency,  two-photon  detuning,  and  the  relative  phase  between  pump  and  Stokes  pulses. 


n 


(  —Qe(t)  cos  A 0  ') 
Qe(t)  sin  A(f) 

V  -m  ) 


we  can  rewrite  equation  43  in  the  compact  form 


(44) 


B  =  Q  x  B . 


(45) 


This  is  the  Bloch  equation,  which  describes  a  precession  of  the  Bloch  vector,  B,  about  the  pseudo 
field  vector,  O,  and  allows  clear,  intuitive  interpretation  of  qubit  dynamics. 
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6.  Alternative  Derivation  of  the  Bloch  Equation 

As  we  showed  in  the  previous  section,  the  density  matrix  can  be  decomposed  using  Pauli  matrixes 
as 

e  =  ±(I  +  B-a).  (46) 

It  is  also  easy  to  verify  that 

B  =  Tr[<f  g]  =  (' Tr[crxg\ ,  Tr[eryo\,  Tr[erzg\ )  =  (u,  v,  w) .  (47) 

Using  the  Pauli  matrix  decomposition  procedure  Tr[aH],  we  can  rewrite  the  Hamiltonian  in 
equation  25  in  the  form 

h  =  |  (n  ■  a) ,  (48) 

where 

f2  =  (—Qe(t)  cos  A(f),  Qe(t)  sin  A(f>,  —S(t))  .  (49) 

Therefore,  the  equation  of  motion  for  the  density  matrix 

q  =  -1-(Hq-qH)  ,  (50) 

takes  the  form 

Q  =  •  Q-  Q-  (0-<t)]  .  (51) 

Using  equation  46,  we  have 

q  =  -l-[(n-*)-(i  +  B-*)^(i  +  B-a)-(n-a)\ 

=  ~[(Q-a)-(B-a)-{B-a)-(Q-a)}.  (52) 

We  know  that  for  the  two  vectors  M  and  N  commuting  with  rr  (see,  for  example,  equation  16.59 
in  reference  (24)) 

(M  ■  B)  ■  (N  ■  (t)  =  M  ■  N  +  ia  ■  (M  x  N)  .  (53) 
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Therefore, 


e  =  ~(Hq-qH) 

=  [(^  •  a)  '  (B  '  *)  ~  (B  ■  a)  ■  (ft  •  <x)] 

=  — i2<r  •  (f2  x  .B) 

4  v  7 

=  (54) 

and  using  equation  47,  we  obtain  the  Bloch  equation 

B  =  Tr[ag\  =  -Tr[a  ■  a  ■  (fl  x  B)]  =  ft  x  B .  (55) 


Equation  55  can  be  written  in  the  matrix  form  as 

B  =  LB , 

where 

/  0  5(f)  £le(f)  sin  A0  ^ 

-E  =  —  5(f)  0  fte(t)  cos  A<j> 

\  —  f2e(f)  sin  A<p  — ft e(t)  cos  Acj)  0  J 


(56) 


(57) 


In  the  Bloch  picture,  the  transformation  of  equation  27  is  the  rotation  of  the  Bloch  vector,  so  that 
the  new  Bloch  vector  is 


B  = 


u  \ 

V 

=  KB  = 

1  „r.  1 

\w  J 


That  gives 

where  ft  =  (fte(t),  0,  —5(f)), 


sin  Acp  cos  Acj)  0 

0  0  1/ 


B  =  ft  x  B  =  LB 


l  u  \ 


w 


L  =  TZLTZ-1 


(  0  5(f)  0  \ 

-5(f)  0  fte(t)  , 

V  0  -fie(f)  0 


(58) 


(59) 


(60) 
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and 


(  cosA</>  sin  A(f>  0  \ 
77.” 1  =  —sin  A  4>  cos  A  <f>  0 

\0  0  1/ 


(61) 


The  transformation  R(t),  equation  30,  in  the  Bloch  picture  is  the  rotation  of  the  Bloch  vector  such 
that 

_  /  \ 

0  1  0 


(  U  \ 

l 

B  = 

V 

=  TZqB  = 

[w  ) 

V 

\  sin [20(f)]  0  cos  [20(f)]  ) 


\w  ) 


where  tan [20(f)]  =  Qe(t)/S(t). 

It  results  in 


B  —  LB  —  TZqTZq  lB  , 


where 


L  —  RqLTZq  i  — 


/  0  A(f)  0  \ 

— A(f)  0  0 

\  o  oo/ 


and 


n 


-l 


C  cos  [20(f)]  0  sin  [20(f)]  ') 

0  1  0 
y  —  sin  [20(f)]  0  cos  [20(f)]  ) 


TZeRe  =  20  (f) 


/  0  0  1  \ 
0  0  0 
V-1  0  0/ 


0(f)  = 


Qe(t)S(t)  -  Qe(t)S(t) 


Neglecting  nonadiabatic  coupling  term  'R.q'JZ0  1  in  equation  63,  we  obtain 

B  =  Cl  x  B , 


(62) 


(63) 


(64) 


(65) 


(66) 


(67) 


(68) 


where  f2  =  (0, 0,  A(f)). 

Note  that  equation  63  can  be  obtained  in  a  different  way  starting  from  equation  31,  constructing 
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equation  for  the  density  matrix,  and  using  7V[er,  g  . 


7.  Evolution  Operator  of  the  Bloch  Vector 


Since  we  already  know  the  exact  form  of  the  evolution  operator  U it)  in  equation  34  and 

|tt(t)>  =  C/(t)|tt(0)>,  (69) 

we  can  easily  construct  the  evolution  operator  for  the  Bloch  vector.  Using  the  definition  of  the 
density  matrix,  we  have 

g(t)  =  |  *(«))<*(«)  |  =  U(i)|4'(0))(4-(0)|!7t(«)  =  U(t)e(0)U'(t),  (70) 

where 


is  the  initial  condition  and 

=  /  cos  (S(t)/2)  —ieiA^  sin  (S(t)/2)  \ 

y  —  ie~lA<t>sin  (S(t)/2)  cos  (<S(f)/2)  J 
=  cos  (S(t)/2)  I  —  i  sin  (S(t)/ 2)  ( elA<t>cr+  +  e~lA<t>crJ) 

=  cos  (S(t)/ 2)  I  +  i  sin  (S(t)/ 2)  (n  •  tr)  =  eiS^n  (T>2  .  (72) 


1  +  Wo  uq  iv  o 
u0  +  iv o  1  -  w0 


(71) 


<?(0)  =  |4r(0))(4r(0)|  = 


£?oo(0)  f?oi (0)  j  _  1 

0io(O)  0n(O)  J  2 


Therefore, 


e(f)  =  t/(f)e(0)(7t(f) 

/  0oo(O)  0oi(O)  \  / 
\  0io(O)  0n(O)  J  \ 


f  cos(S'(f)/2)  ielA<t>  sin  (S(t)/ 2)  \ 

y  ie_*A^sin  (S(t)/2)  cos  (<S(f)/2)  J 

cos(S'(f)/2)  —  ielA^sin  (S(t)/2)  \ 

— ie~lA<t>sm  (S(t)/2)  cos  (<S(f)/2)  J 


(73) 
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or,  exactly, 


Qoo(t)  =  Qoo(0)  cos2  (S(t) /2)  +  Qn(0)  sin2  (S(t)/2) 

+  ^  sin  (£(£))  (gw(0)elA4,  -  ^oi(O)e_jA0)  ,  (74a) 

goi(t)  =  ~etA0  sin  (S^i))  (goo(0)  -  £n(0))  +  g10(0)e2lM  sin2  (S(t)/2) 

+  0oi(O)  cos2  (S'(t)/2)  ,  (74b) 

gw(t)  =  ^e_*A<^  sin  (S(t))  (^oo(0)  -  £n(0))  4-  Qm(^)e~2lA4>  sin2  (S(t)/2) 

+  0io(O)  cos2  (S(t)/2):  ,  (74c) 

Qn(t)  =  0oo(O)  sin2  (S(t)/2)  +  pn(0)  cos2  (S(t)/2) 

+  ^  sin  (£(£))  (^oi(O)e"*A0  -  pio(O)e*A0)  .  (74d) 


Taking  into  account  equation  41,  we  obtain  for  the  Bloch  vector  components 

u(t )  =  sin  (S(t))  sin  (A0)  wq  +  (cos2  (S(t)/2)  +  sin2  (S(t)/2)  cos  (2A0))  uq 

—  sin2  (S(t)/2)  sin  (2A(f>)  v0  ,  (75a) 

v(t)  =  sin  (S(t))  cos  (A0)  wq  —  sin2  (S(t)/2)  sin  (2A0)  uq 

+  (cos2  (S(t)/2)  —  sin2  (S(t)/2)  cos  (2A0))  v0  (75b) 

w(t)  =  (cos2  (S(t)/2)  —  sin2  (S(t)/2))  w0 

—  sin  (S'(f))  («o  sin  (A0)  +  n0  cos  (A0))  ,  (75c) 

which  can  be  presented  in  the  matrix  form  as 

/  C2  +  S2  cos  (2A0)  — 52sin(2A0)  2C-5sin(A0) 

B(t)  —  I  —  52sin(2A0)  C2  —  S2  cos  (2A0)  2C-5cos(A0)  -B(O),  (76) 

\  —  2C-5sin(A0)  — 2C-5cos(A0)  C2  —  S2  ) 

where  C  =  cos  (S(t) /2),  S  =  sin  (S(t)/2). 
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8.  Ultrafast  Qubit  Rotations  Using  Geometrical  Phase 


At  this  point,  we  are  ready  to  discuss  implementation  of  the  single  qubit  gates  since  we  have 
obtained  the  analytic  solution  for  the  qubit  wave  function  and  constructed  the  evolution  operator 
in  the  Bloch  vector  representation.  A  universal  set  of  quantum  gates  has  been  intensively 
discussed  in  the  literature  related  to  the  universality  in  quantum  computation  (9-72).  To  perform 
quantum  computation,  we  must  have  two  major  building  blocks  at  our  disposal:  arbitrary  unitary 
operations  on  a  single  qubit  and  a  controlled-NOT  operation  on  two  qubits.  Here  we  address 
only  single  qubit  manipulation. 

To  demonstrate  arbitrary  geometric  operations  on  a  single  qubit,  we  use  the  Bloch  vector 
representation  discussed  in  the  previous  section.  Since  any  unitary  rotation  of  the  Bloch  vector 
can  be  decomposed  as  (9,  13) 


U  =  eUxoRz(a1)Ry(a2)Rz(a3) ,  (77) 

where  R,  =  eia(Ti  (i  =  y,  z)  are  the  rotation  operators,  we  need  to  demonstrate  rotations  of  the 
qubit  Bloch  vector  about  the  z  and  y  axes  by  applying  various  sequences  of  external  pulses.  The 
decomposition  in  equation  77  plays  an  important  role  in  circuit-based  quantum  computing,  as  it 
shows  explicitly  that  two  single-qubit  operations  allow  us  to  create  the  arbitrary  state  of  the  qubit. 
Here  we  show  how  this  can  be  accomplished  by  controlling  the  parameters  of  the  external  pulses, 
which  are  defined  by  the  explicit  form  of  the  evolution  operator  (see  equation  34).  There  are  two 
distinct  ways  of  the  implementation  depending  on  which  part  of  the  total  qubit  phase  we  employ: 
dynamical  or  geometrical  (14).  Quantum  gates  relying  on  geometrical  quantum  phases  are  called 
holonomic  gates  and  they  are  expected  to  be  robust  with  respect  to  noise  (15,  18,  19). 

To  implement  the  rotation  of  the  Bloch  vector  about  z  axis  (the  phase  gate)  based  on  the 
geometrical  phase,  we  can  use  the  evolution  operator  of  the  resonant  qubit,  equation  34.  The 
product  of  two  evolution  operators  corresponding  to  the  sequence  of  two  n  pulses  with  the 
relative  phase  A<p  =  ip  +  vr,  gives 

(  ei<f>  0  \ 

Rz(ip)  =  0  =  .  ,  (78) 

\0  e  I 

where  the  first  subindex  of  U  indicates  the  pulse  area,  S(T),  and  the  second  one  indicates  the 
phase,  A(f>. 
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Figure  5.  The  Bloch  vector  trajectory  for  the  qubit  state  |0)  in  panel  (a)  and  the  qubit  state  1 1) 
in  panel  (b)  generated  by  the  sequence  of  two  n -pulses  with  the  relative  phase  ip  +  vr. 


Figure  5  shows  the  Bloch  vector  trajectories  of  the  qubit  basis  states  |0)  and  |1),  which 
correspond  to  the  angles  (3  =  0  and  (3  =  tx  in  equation  41  and  the  Bloch  vector  initially  pointing 
in  z  and  —2  directions  while  the  vector  f2i  =  (— De,  0,  0)  is  pointing  in  —  x  direction.  For 
simplicity  we  chose  A0  =  0  for  the  first  n -pulse.  The  fist  7r-pulse  flips  the  population  to  the  state 
|1)  (|0));  correspondingly,  the  Bloch  vector  turns  about  the  effective  field  vector  Oi  (about  the  x 
axis),  and  it  stays  in  the  y,  z  plane  all  the  time  and  points  in  the  —  z  (z)  direction  at  the  end  of  the 
pulse.  Due  to  the  second  7r -pulse,  the  population  is  transferred  back  to  the  initial  state  |0)  (|  1)); 
therefore,  the  Bloch  vector  returns  to  its  original  position  pointing  along  the  0  (—z)  axis. 

However,  since  we  chose  A0  =  ip  +  7r  for  the  second  7r-pulse,  the  pseudo-field  vector  is  rotated 
counterclockwise  by  the  angle  <p  +  ix  in  the  x,  y  plane,  02  =  (fie  cos  ip,  —Qe  sin  ip,  0),  and  the 
Bloch  vector  moves  in  the  plane  perpendicular  to  the  x,  y  plane  and  has  the  angle  7t/2  —  a 
(— 7t/2  —  a)  with  the  x,  z  plane. 

The  Bloch  vectors  representing  a  pair  of  orthogonal  basis  states  |0)  and  |1)  follow  a  path 
enclosing  correspondingly  solid  angles  of  2<p  and  —  2p.  The  geometrical  phase  is  equal  to  one 
half  of  the  solid  angle,  which  means  the  basis  states  |0)  and  |1)  gain  phases  <p  and  —p  and  the 
evolution  operator  takes  the  form  of  the  phase  gate,  equation  78,  with  the  relative  phase 
controlling  the  phase  of  the  gate. 

The  rotation  operator  about  the  y  axis  can  be  constructed  using  three  pulses.  The  first  and  third 
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pulse  is  7r/2-pulse  with  A 0  =  0,  while  the  second  pulse  is  tt -pulse  with  the  relative  phase  tt  +  <p. 
It  is  easy  to  show,  using  equation  34,  that  this  three-pulse  sequence  results  in 

Ry(ip)  =  U^oU^+yU*  ;0  =  e^.  (79) 


Figure  6.  The  Bloch  vector  trajectory  for  the  qubit  state  | i)  in  panel  (a)  and  the  qubit  state  |  —  i) 
in  panel  (b)  generated  by  the  sequence  of  two  tt / 2-pulses  and  one  tt  pulse  with  the 
relative  phase  tp  +  tt. 


To  demonstrate  the  geometrical  nature  of  the  Ry{<p)  operation,  we  use  the  fact  that  it  creates  the 
relative  phase  between  the  qubit  basis  states  |  ±i)  —  (|0)  ±  i\l))/\/2  (see  Appendix  for  the 
eigenvectors  of  the  Pauli  matrices).  In  the  Bloch  representation,  these  states  have  the  form 

I  ±  i)  =  cos  (0  |0)  +  e±l 5  sin  (0  |1) ,  (80) 

which  are  two  vectors  defined  by  the  angles  (5  =  tt / 2  and  a  =  ±7t/2  and  pointing  in  the  y  and  —y 
directions,  as  shown  in  figure  6.  The  trajectory  of  the  Bloch  vector  representing  the  states  |  ±  i) 
is  shown  in  figure  6.  The  pseudo-field  vectors  f2i  and  03  are  defined  by  the  effective  Rabi 
frequencies  of  the  first  and  third  pulses  and  are  pointing  in  the  —x  direction  since  A©  =  0.  The 
second  pseudo-field  vector  02  is  rotated  counterclockwise  by  the  angle  <p  +  tt  in  the  x,  y  plane 
same  as  in  the  case  above.  The  initial  Bloch  vector  is  pointing  in  the  y  (—y)  direction.  The  first 
7r/2-pulse  rotates  the  Bloch  vector  about  Oi  to  the  position  of  the  state  |1)  (|0)).  The  second 
pulse  flips  the  direction  of  the  Bloch  vector.  The  third  7r/2-pulse  returns  the  Bloch  vector  to  its 
original  position.  The  Bloch  vector  and  the  pseudo-field  vector  are  orthogonal  during  the  whole 
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evolution.  Similar  to  the  previous  case,  we  observe  that  the  basis  states  |«)  and  |  —  i)  follow  a 
path  enclosing  correspondingly  solid  angles  of  2p  and  —2 p.  Therefore,  they  gain  the  relative 
phase  2p,  which  is  the  geometrical  phase  defined  by  the  relative  phase  between  pulses.  It  is  easy 
to  show  that  the  phase  gate  in  the  |  ±  i)  basis  is  equivalent  to  the  Ry(p)  gate  in  the  |0),  |1)  basis. 


9.  Rotation  in  the  Bloch  Representation 


Rotation  operations  in  the  Schrodinger  picture  are 


I  tt>  =  Ri|*o>, 


where  i  =  {x,  y,  z},  |T0)  is  the  initial  wave  function, 

R  =  ei(pcrx  = 


1  =  g-ilpo*  = 


Ry  =  eitpcTy  = 


R-l  =  e-VP<Ty  = 


cos  ip 

i  sin  p  \ 

i  sin  (p 

cos  p  J 

cos  p 

—i  sin  <y 

— i  sin  p 

cos  p 

cos  p 

sin  p  \ 

—  sin  p 

cos  p  J 

1  cos  ip 

—  sin  p 

l  sinyj 

cos  p 

Ry  =  e ^  = 


evp  0 
0  e~vp 


1  =  e_Up<rx  = 


e~ltp  0 
0  evp 


(81) 


(82) 


(83) 


(84) 


(85) 


(86) 


(87) 
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Therefore  for  the  density  matrix,  we  have 


Q  =  |*)<*|  =  Ri\V  o)(^o|Kr1  =  RiOoRi1 ,  (88) 

where  go  =  To)  (To  is  the  initial  density  matrix. 

Using  equation 


B  =  Tr[ag]  =  ( Tr[axg],Tr[ayg ],  Tr[crzg ])  =  (u,  v,  w ) ,  (89) 

we  find  the  following  expressions  for  the  transformation  of  the  Bloch  vector  components 

u  =  Tr[crxRg0Rr x] ,  (90) 

v  =  Tr{<jyRg0R-x} ,  (91) 

w  =  Tr[crzRg0Rr1] .  (92) 


Substituting  equations  82  -  87  into  equations  90  -  92,  we  obtain 


/I 


TZr  = 


0 


0  \ 


0  cos  (2 Lp)  sin  (2 ip) 

\  0  —  sin  (2 tp)  cos  ( 2tp )  ) 


IZy  ~ 


(  cos  (2 ip)  0  —  sin  (2 ip)  ^ 
0  1  0 

\  sin  (2 ip)  0  cos  (2 ip)  ) 


TZ?  = 


(  cos  (2(p)  sin  (2(f)  0  ') 

—  sin  (2(f)  cos  (2(f)  0 

0  0  1/ 


(93) 


(94) 


(95) 


From  equation  76,  we  see  that  the  evolution  operator  in  the  Bloch  representation  is 


m 


(  C2  +  S2  cos  (2A cj)) 
—S2  sin  (2Ao) 

\  -2C-S  sin  (A0) 


— S2  sin  (2A0) 
C2  —  S2  cos  (2A0) 
—2  C  ■  S  cos  (A0) 


2  C  ■  S  sin  (A0)  ^ 
2  C  ■  S  cos  (A 4>) 
C2-S 2  j 


(96) 
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Therefore,  to  demonstrate,  for  example,  the  Rz  rotation  in  the  Bloch  picture,  we  see  that  the  first 
7r  pulse  with  the  relative  phase  A <f>  =  0  gives  for  the  evolution  operator 


ff- ;  | ) 


/  1  0 
0  -1 
V  o  o 


0  \ 
0 

-1/ 


(97) 


The  second  n  pulse  with  the  relative  phase  A©  =  tt  +  ip  gives  for  the  evolution  operator 


U 


7T;7T  -\-<p 


(  COS  (2(7T  +  p)) 
-  sin  (2(7t  +  p)) 

V  0 


—  sin  (2(7t  +  </?))  0  ^ 

—  COS  (2(7T  +  93))  0 

0  -1  / 


Finally,  the  sequence  of  U7r-n+ip  and  £4-;o  results  in 


(98) 


( 

COS  (2(7T  +  p)) 

-  sin  (2(7 r  +  p)) 

0  \ 

< 

1 

0 

0  \ 

^7r;7r+(/?^7r;0 

—  sin  (2(7t  +  p)) 

—  cos  (2(7 r  +  p)) 

0 

0  -1 

0 

V 

0 

0 

-1 ) 

\ 

0 

0 

-1/ 

C  cos  (2y))  sin  (293)  0  ^ 

—  sin  (2p)  cos  (2p)  0 

V  0  01/ 


(99) 


which,  of  cause,  coincides  with  equation  95.  Similarly,  one  could  consider  the  Ry  rotation  in  the 
Bloch  picture. 


10.  Generalization  of  the  Single-Qubit  Operation  Using  Bright-Dark  Basis 


In  the  previous  sections,  we  have  considered  several  excitation  schemes  of  the  three-level  system 
and  discussed  a  possible  implementation  of  single-qubit  gates.  It  was  shown  that  all  possible 
qubit  states  can  be  created  in  a  controllable  fashion  using  a  couple  of  completely  overlapped  laser 
pulses,  fip(t)  =  (15(f).  In  this  section,  we  present  a  bit  more  general  solution,  which  allow  some 
additional  flexibility  in  terms  of  the  ratio  of  the  pump  and  Stokes  pulse  amplitudes.  Again,  we 
address  here  the  coherent  Raman  excitation  in  a  three-level  A-type  system  consisting  of  the  two 
lowest  states  of  electron  spin  |0)  =  |  —  X)  and  |0)  =  A")  coupled  through  an  intermediate  trion 
state  | T)  (see  figure  1)  and  assume  that  the  trion  state  is  far  off-resonance  with  the  external  fields. 
In  addition,  we  restrict  our  consideration  to  the  non-impulsive  regime,  and  can  then  put 
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nP+(t)  =  ns+(t)  =  0  in  the  equation  6  so  that  the  Hamiltonian  takes  the  form 

(100) 

(101) 

(102) 

where  Ap  —  u op  —  up,  Au  =  up  —  us,  fip+  =  (2p0(f)e^pO),  Qs+  —  fls0(f)e^sO)5  anc[ 

4>P,s(t )  =  0p,S  +  oip^st2  /  2. 


H  —  U rwaHUrwa  iftURWAUpcwA 


/ 

0 

0 

Hp+ 

(  ° 

0 

0  \ 

0 

—2c op 

+  h 

0 

—Au 

0 

V 

0* 

“p+ 

O* 

“5+ 

— 2c up  y 

\o 

0 

—Up  J 

/  0  0  0P+  \ 
0  2(Acu  —  cup)  (I5+ 

V  fip+  ^s+  -2Ap  / 


Let  us  consider  the  case  when  fipo(t)  =  H0(f)  cos  i?  and  fiso(f)  =  O0(f)  sin-#;  the  same 
time-dependent  envelope  for  the  pump  and  Stokes  Rabi  frequencies  while  the  mixing  angle  -# 
controls  the  ration  between  the  maximum  of  the  Rabi  frequencies.  In  this  case,  we  make 


transformation  to  the  bright-dark  basis  (  'k)  =  jap(t),  ap(f),  61 

;<)} 

=  7£M|(k))  using  the 

following  matrix 

(  e  l(tp  cos ‘t)  e  l(t’s  sin  i) 

0  \ 

= 

sjn  {j  ei<f>p  cos  $ 

0 

(103) 

^  0  0 

1 

/  el^p  cos  t)  —  sin  # 

0  3 

7?-1  — 
'Hd  — 

ei<t>s  sjn  ■{}  e-i<t>p  cos  $ 

0 

(104) 

\  0  0 

ly 

In  the  bright-dark  basis,  the  Hamiltonian  takes  the  form 

h  =  nhdHK,;} 

/  2  SS2  2 5SCe-^+  -f+(t)n0(t)  \ 

=  |  2  SSCe^  2  SC2  e**+£-(t)SCn0(t)  , 

V  -C+(t)n0 (t)  e-*+£(t)SCno(t)  2AP 

where  S  =  sin  t),  C  =  cos  t),  4>+  =  4>p  +  4>s ,  £+(£)  =  e*"‘2//2C2  +  etf3t2/2S2  , 
f_(t)  =  eiat2/2  -  e^2/2. 


(105) 
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In  the  resonant  case,  5  =  0,  we  obtain 


H 


lZbdHlZhd 

h  (  0  0  -Z+{t)too(t)  \ 

^  0  0  e^+^(t)scn0(t) 

V  -e+(t)n0(t)  e~^+C(t)SCn0(t)  2  A p  ) 


(106) 


In  the  case  of  equal  chirp  rates,  a  =  (3,  we  have 


2 SSCe-V*  -eiat2/H \(t)  \ 

2  SC2  0 

0  2A  P 

In  the  two-photon  resonance  case,  5  =  0,  and  equal  chirp  rates,  a  =  /3,  we  have 

H  =  nbdHTZ^ 

0  0  -eiat2/2tt0  (t)  \ 

0  0  0 
-e~iat2/2n0  (t)  0  2  A  P 


H  =  nbdHU^ 

(  2  3S2 

h 

=  -  2  5SCel*+ 

V  -e~iat2/H \{t) 


(107) 


(108) 


We  can  see  from  equation  108  that  a  nit)  =  aB(  0)  and  we  are  left  with  the  system  of  two 
differential  equations  for  the  probability  amplitudes  aB{t)  and  b(t).  Making  the  adiabatic 
elimination  of  the  excited  state  | T)  (assuming  that  b(t )  ~  0),  we  have 

bit)  =  A-n0(t)e-iat2/2aB(t),  (109) 

and 

iaB(t)  =  —  (110) 


The  solution  of  the  equation  1 10  is 


aB(t) 


aB(0)e~4 


■  fo  Q ' 


(111) 


Therefore,  the  evolution  operator  for  the  two-photon  resonant  excitation  by  the  equally  chirped 
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pulses  under  condition  of  the  adiabatic  elimination  of  the  trion  state  has  the  following  form 


/  g4Ap  /o  Q  Q  \ 


u(t)  =  n, 


-1 

bd 
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o 
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1  o 
o  1 


n, 


bd 
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=  |  el(t>s  sin d  e~l<^p  cost!)  0 

0  0  1/ 
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0  1  0 
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o 


/  e'lS(t)  cos2  d  +  sin2  d  e*A</>(e?5(t)  _  i)  sinr?cosr?  0\ 


iA(f)  ^ 


1)  sin  d  cos 'd  cos2  d  +  e lS^  sin2  d 


V 


0  0 
where  S(t)  =  J*  Q%(t')dt' /(4AP),  and  A 0  =  0p  —  0s. 

For  the  qubit  states,  the  evolution  operator  can  be  written  as 


1  ) 

(112) 


1/ 


U(t)  =  e 


-  JS(t)/ 2 


elSdy2 cos2 ’d  +  e  tS^/2sin 2  d  eiA<t>^eiS(t)/2  _  e  iS{t)/2^  sjn$cos$ 
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iS(t)/2  (  cos(5'(f)/2)  +  i  sin  (5(f) /2)  cos  (2d)  ietA^sin  (S(t)/2)  sin  (2d) 

y  ie~tA<f>sin  (S(t)/2)  sin  (2d)  cos  (S(t) /2)  —  i  sin  (S(t) /2)  cos  (2d) 
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(113) 


_  eiS(t)/2e-iS(t)n-cr/2 


where 


n  =  (—  cos  (A0)  sin  (2d) ,  sin  (A0)  sin  (2d) ,  —  cos  (2d)) 


(114) 


11.  Electron  Spin  in  a  Quantum  Dot  as  a  Qubit 


In  the  previous  sections,  we  developed  several  methods  of  an  arbitrary  manipulation  of  a  qubit 
wave  function  using  the  geometric  phase.  Now  we  apply  the  proposed  scheme  to  electron  spin 
states  in  a  charged  quantum  dot.  Due  to  quantum  confinement,  the  state  of  the  electron  can  be 
expressed  as  a  product  of  the  Bloch  function  and  an  envelope  function,  which  has  a  typical  scale 
of  the  quantum  dot  size,  a  few  nanometers.  The  energy  level  structure  and  optical  selection  rules 
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have  been  discussed  in  the  literature  (1-4).  A  commonly  accepted  energy  level  structure  is 
comprised  of  four  levels,  the  two  electron  spin  states  and  two  trion  spin  states.  Figure  7  shows 
two  arrangements  of  the  energy  levels  and  polarization  selection  rules,  which  provide  a  possibility 
of  optical  control  for  the  electron  spin  qubit. 


Figure  7.  Optical  selection  rules  in  different  bases. 

The  usual  experiment  of  the  electron  spin  control  is  performed  at  low  temperature  (~  1  K).  An 
external  magnetic  field  in  the  Voigt  configuration  (of  order  2  -F  7  T)  is  applied  along  the  x-axis, 
perpendicular  to  the  sample  growth  direction,  the  z  axis.  Zeeman  splitting  of  the  electron  and 
trion  spin  states  are  on  the  order  a ;e  =  10  meV  and  ay  =  10/xeV,  correspondingly.  At  these 
conditions,  taking  into  account  the  optical  selection  rules,  the  four-level  system  can  be  considered 
as  a  double  A  system  (1-4,  25).  This  coupling  scheme  is  shown  in  figure  7(b),  where  we 
indicated  by  “77”  and  “V”  the  optical  field  couplings  with  the  orthogonal  polarization.  The 
shown  coupling  scheme  is  in  the  so-called  “x  basis.” 

An  alternative  arrangement  is  depicted  in  figure  7(a).  In  this  case  the  mixed  basis  is  used,  where 
the  electron  spin  states  are  in  the  x  basis  while  the  trion  states  are  in  the  z  basis  (1—4,  25).  Using 
the  two  a+  or  o~  polarized  fields,  one  can  couple  the  electron  spin  states  V_)  and  X+)  as 
shown  in  the  figure  7(a).  This  is  the  case  where  our  three-level  model  can  be  implemented.  The 
corresponding  Hamiltonian  has  the  form  obtained  in  equation  16.  Assuming  large  detunings  of 
the  pump  and  Stokes  field  frequencies  from  the  transition  frequencies  to  the  trion  state  T),  after 
the  adiabatic  elimination  of  the  trion  state,  for  the  case  of  completely  overlapped,  identically 
linearly  chirped  pump  and  Stokes  pulses,  we  obtain  the  following  Hamiltonian 


(115) 
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where  5  =  ue  —  Au,  and  Au  =  uP  —  us,  f2e(i)  =  QP0(t)QS0(t)[  1  +  e_^A<?i+AaJ^]2/(2A)  is  the 
effective  two-photon  Rabi  frequency,  A  =  ut  —  top,  and  A<j)  =  Op  —  Os-  Note,  that  here  we 
have  used  /x0t  ~  /Jit- 


t/r0 


Figure  8.  The  population  dynamics  of  the  resonant  qubit  states  with  (dotted  lines) 
and  without  (solid  lines)  adiabatic  elimination  of  the  trion  state  in  the 
three-level  system.  The  excitation  is  generated  by  the  sequence  of  two 
pairs  of  7r -pulses  (Gaussian  pulse  envelops)  with  the  relative  phase 
Aip  =  7t/2. 


We  observe  almost  perfect  resemblance  between  equation  115  and  equation  25.  The  only 
difference  is  that  the  effective  two-photon  Rabi  frequency  contains  the  oscillating  term  with  a 
frequency  Au;.  Now,  we  show  that  the  contribution  of  this  oscillating  term  can  be  neglected 
when  the  pulse  duration  is  longer  than  Au;-1  ~  u~l.  To  demonstrate  this  and  justify  the 
procedure  of  the  adiabatic  elimination  of  the  trion  state,  we  numerically  solve  the  time-dependent 
Schrodinger  equation  with  the  Hamiltonian  in  equation  16.  That  is,  we  compare  our  analytic 
solution  with  the  exact  solution  of  the  Schrodinger  equation  without  the  adiabatic  elimination 
approximation.  An  example  of  the  comparison  is  shown  in  figure  8. 

The  example  in  figure  8  shows  the  dynamics  of  the  population  of  the  resonant  qubit  excited  by  a 
sequence  of  two  pairs  of  linearly  chirped  pulses  with  the  relative  phase  between  the  pairs  being 
< j)  =  7t/2.  The  parameters  of  the  excitations  are  the  transform-limited  pulse  duration  To  =  100  fs, 
the  pulse  area  of  the  pump  and  Stokes  pulses  is  equal  to  7 r,  liner  chirp  rate  a' / r02  =  20,  Ar0  =  6, 
ujeT o  =  1,  the  ratio  between  the  maximum  Rabi  frequency  and  the  single-photon  detuning 
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f20 /A  =  0.77.  We  observe  a  negligibly  small  amount  of  population  in  the  trion  state  | T)  at  the 
intermediate  time.  The  presented  comparison  also  demonstrates  a  reasonable  agreement  between 
the  proposed  control  schemes  and  the  exact  numerical  solution.  Note  that  the  total  time  for  the 
qubit  operation  in  figure  8  is  on  the  order  of  25  ps,  which  is  much  shorter  than  the  typical  life  time 
of  the  trion  state  ( 1 ,  2,  25)  as  well  as  the  time  scale  of  other  forms  of  decoherence,  such  as  that 
induced  by  the  electron-phonon  interaction. 


12.  Conclusion 


We  presented  the  analytic  expression  of  the  evolution  operator  of  the  electron  spin  in  a  quantum 
dot,  which  provides  a  clear  geometrical  interpretation  of  qubit  dynamics.  Using  the  analytic  form 
of  the  evolution  operator,  we  proposed  a  set  of  single-qubit  rotations  that  are  solely  based  on  the 
geometrical  phase.  Our  proposal  combines  the  pulse  area  control  with  the  adiabaticity  by  using 
chirped  pulses.  To  estimate  the  time  scale  of  the  proposed  operation,  we  can  use  the  100  fs 
pulses,  which  become  picosecond  pulses  after  linearly  chirping.  That  amount  of  chirping  is 
sufficient  to  provide  adiabatic  excitation  (6,  7,  20)  and  can  be  readily  produced  experimentally 
using  commercially  available  laser  systems.  Using  parameters  of  the  dipole  moments 
/j,ot,ti  ~  200  D  available  in  InGaN/GaN  (27),  GaN/AIN  (22,  23)  quantum  dots,  and  detuning 
A  =  5  meV,  we  estimate  the  peak  amplitude  of  the  pulses  on  the  order  of  106  —  107  V/m.  Note 
that  the  demonstrated  adiabatic  manipulation  of  a  qubit  using  only  the  geometric  phase  has  some 
advantages,  since  it  reduces  the  requirements  of  perfect  tuning  of  the  control  field  parameters  and 
is  significantly  more  robust  against  noise  (75,  18,  19). 
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Appendix.  Eigenvectors  of  the  Pauli  matrices 


The  eigenstates  of  crx  are 


I  j)  = 


+  |1)) 
V2 


=  cos  Q)  |0)  +sin  Q)  |1)  = 


1  /  1 


V2  \  1 


(A-l) 


(A-2) 


Correspondingly, 


0  and  (3  —  — ,  a  —  0, 


0  and  (3  =  — ,  a  =  n. 

0  / 


(A-3) 


(A-4) 


The  eigenstates  of  cry  are 
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(A-6) 


Correspondingly, 


|?;)  <-»•  B=  1  and  p=^,a  =  ^, 


—  i)  B  = 


.  n  n  3vr 

and  f3  —  —,a  —  — . 
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The  eigenstates  of  <rz  are 


|0>  = 
|1>  = 


(A-9) 
(A- 10) 


Correspondingly, 


|0> 

|1) 
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(  0  \ 
o 


and  (3 


0,  a  is  arbitrary, 


V  1  / 

(  0  \ 

o 

V  -i  / 


and  (3 


n.  o.  is  arbitrary. 
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